Does the melting behavior in ultrathin metalhc nanowires begin from the surface ? 
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Surface melting temperature is well-known significantly 
lower than the bulk melting point. But we find that the inte- 
rior melting temperature in ultrathin nanowires is lower than 
that of the surface melting . The thermal stability of heli- 
cal multi-walled cylindrical gold nanowires is studied using 
molecular dynamics simulations. The melting temperature 
of gold nanowires is lower than the bulk value, but higher 
than that of gold nanoclusters. An interesting interior melt- 
ing is revealed in the gold nanowires and the thermodynamics 
is closely related to interior structures. The melting starts 
from the interior atoms, while the surface melting occurs at 
relatively higher temperature. We propose that the surface 
melting represents the overall melting in ultrathin metallic 
nanowires. 
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The melting behavior of nanoparticles and nanorods 
have been demonstrated dramatically different jiom the 
bulk both experimentally and theoreticallyuilj. The 
melting process of a crystalline starts from the surface 
layer and propagates into the interior. Thus, the surface 
melting temperature is significantly lower than the bulk 
melting point. Similarly, people may ask whether the sur- 
face melting temperature is lower than the overall melt- 
ing temperature in clusters and nanowires. Berry consid- 
ers that "dynamic coexistence" or surface melting hap- 
pens in the iiielting process of small clusters before the 
total meltingp. For crystalline nanowires, Tosatti found 
that surface melting temperature for Pb |Wires was also 
lower than the total melting temperaturecl. Experimen- 
tally, surface melting is involved Jn the melting process 
of nanoparticles and nanorodsElLl. Especially, Schmidt 
found a broad peak for heat capacity of Najjg clusters, 
maybe implying the coexistence of |Solid-like and liquid- 
like phases before the total meltinga. Two major effects 
are responsible for these different melting behavior in 
nanoparticles and nanorods. One is the large surface-to- 
volume ratio in these nanostructures. The other is quan- 



tum confinement effect in finite size systems. Surface 
atoms have fewer nearest neighbors and weaker binding, 
which may induce an earlier surface melting behavior. 
On the other hand, a close relationship between the pfilt ,- 
ing and the structure features is found in the clustertillla. 

Recent studies demonstrated that ultrathin metal- 
lic nanowires have quite different structural and prop- 
erties from r-those of bulk, clusters and crystalline 
nanowiregl3~E3. The helical multi-walled cylindri- 
cal structures have been widely found in l-3nm size 



range of meta' 
theoreticallyli3^ 
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lowires both experimentally and 
This kind of novel structure will 
bring some bizarre melting features different from the 
above-mentioned systems. To our knowledge, less efforts 
are focused on their thermodynamics so far, although this 
kind of uJlpithin metallic nanowires have attracted great 
interestsEJ'B. 

Furthermore, the ultrathin nanowire has some charac- 
teristics similar to the cluster, crystalline wire and bulk. 
It may provide an opportunity for comprehensively un- 
derstanding these types of matter states and their rela- 
tions. For example, the surface-to-volume ratio is a non- 
zero value in these nanostructures, while it approaches 
toward zero in the bulk. Surface and core atoms are ex- 
pected to play different role during the melting process. 
However, it is difficult to definitely distinguish the over- 
all melting and surface melting of clusters experimentally 
because the cluster's signal spectrum will not exist when 
the surface melting occurs. Therefore, we can employ 
the ultrathin gold nanowire as a representation to ex- 
plore these problems. 

In this letter, the thermal stabilities of gold nanowires 
with helical multi-walled cylindrical structures are stud- 
ied by using molecular dynamical (MD) simulations. 
We start from the optimized structures from previous 
worksllj, which were imaged by electron microscopyllZI. 
The interactipa between gold atoms is described by a 
glue potentialEJ and the periodic boundary condition is 
applied along the wire axis direction to model the wire 
with sufficient length wire. Xhc length of supcrcell is 
chosen with the same as Ref.Ej, which is a reasonable 
scheme to attain the helical structures in the nanowire. 
To characterize the thermal behavior of nanowires, we 
monitor the root-mean-square (rms) fluctuation of the 
interatomic bond distances S defined by 
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and the heat capacity C per atom, which is related to 
the energy E fluctuation, by the relation 
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where r(i,j) denotes the distance between the nuclei i 
and j, n is the total number of the atoms in the nanowire, 
ks is the Boltzman constant, (...) indicates the thermal 
statistical averages in the canonical ensemble after equi- 
libration. 

The constant temperature molecular dynamics (MD) 
method of Nosecj is employed to exploit the thermal 
properties of gold nanowires. The MD time step is cho- 
sen as 2.15xl0~^^s. The initial 10^ steps are used to 
bring the system into equilibration and further 10^ steps 
are used to record the thermal average of the physical 
quantities. We study five representative helical struc- 
tures of gold nanowires with up to triple shells: single 
atom centered, double-chain, trigonal, tetragonal and six 
atoms parallelogram centered noted as SI, S2, S3, S4, S6 
in Fig.l. Table I presents the starting melting tempera- 
ture Tini and the overall melting temperature Tm of these 
structures. The overall melting temperature Tm reflects 
the stability of the whole system, while the starting melt- 
ing temperature T^i describes the early stage of melting 
and sensitively depends on wire structures. The over- 
all melting temperature T,„ is estimated from the curves 
of rms bond length fluctuation S, heat capacity C and 
binding energy E as functions of temperature. 




FIG. 1. Morphology of some representative gold nanowires 
with triple shells. 

As shown in Table I, the overall melting temperatures 
r„j for nanowires are almost the same for all the wires 
and are slightly lower than the bulk melting point. Ex- 
perimentally, Liu et al. have rspprted the melting point 
of Pt nanowire is about 400°Clli. Lee et al. have also 



found the melting temperature of 4.6nm Pd nanowire 
is just 300,°C, which is much lower than the bulk value 
(1445°C)ll3. Similar phenomenon is also found in the 
case of metal clusters. The depression of melting point 
can be attributed to the low dimension and large surface- 
to-volumc ratio in these nanostructures. However, the 
melting temperatures obtained for gold nauQwires are 
much higher than those of gold nanoclusterstj This 
may be understood by the tightened helical structures 
in nanowires. On the other hand, the starting melting 
temperatures Ti^i are different for wires with different 
interior structures (see Table I). The wire SI (as de- 
scribed by multi-walled structural index: 18-12-6-1) starts 
to melt at a rather lower temperature than the other 
ones, while the wire S4 (21-15-9-4) has relatively higher 
starting melting temperature. These indicate that the 
different starting temperature Tint for different nanowires 
are related to their different interior structures. Similar 
overall stability for all the wires may come from their 
common multi-walled helical packing, despite different 
interior structure. Moreover, since the Ti„i is related to 
the stability of interior atomss, the rather lower start- 
ing melting temperature T^i implies that the interior 
melting behavior happens during melting process of the 
helical gold nanowires. 

We discuss the structural evolution of gold nanowires 
during the melting process. Fig. 2 gives several snap shots 
taken from the structural trajectories of the wire SI at 
different temperatures. It is interesting to note that the 
interior atoms diffuse along the wire axis direction at a 
rather low temperature. As shown in Fig. 2, the center 
atoms firstly move along the wire at 300K (Fig. 2a). The 
helical structure of the outmost shell is almost invariant. 
With the raise of temperature, the center atoms continue 
moving away from the wire. The atoms in the first shell 
(from interior to outer) thus have fewer nearest neighbors 
and begin to diffuse along axis direction (Fig. 2b). Simi- 
larly, the atoms in the second shell are also involved into 
the migration at higher temperature T — 900K (Fig. 2c). 
However, the helical structure of surface shell still ex- 
ists. When the temperature is high enough, the surface 
atoms in the outmost shell also start to migrate and the 
hehcal structure of surface is broken at HOOK, leading 
to the overall melting (Fig. 2d). Therefore, we conclude 
that interior atoms diffuse prior to surface atoms and no 
surface melting takes place before the overall melting in 
gold nanowires. 

Table I. Melting temperature for different interior struc- 
ture of Au nanowire. 
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FIG. 2. Structure evolution of gold nanowire 81(18-12-6-1) 
with elevated temperature. 

To further illustrate the argument of the interior melt- 
ing and distinguish the role of surface and core atoms 
on the melting behavior, Fig. 3 plots the rms bond length 
fluctuation of the surface atoms (5s), core atoms (Sc) and 
the total atoms {S) of SI wire as functions of tempera- 
ture. Obviously, the rms bond length fluctuation of the 
core atoms have similar trends to that of the whole wire, 
but dramatic different from the surface. In the temper- 
ature range of 350 — lOOOK, the Sg of the surface atoms 
is very small and almost invariable, while 6c and 6 have 
substantial fluctuations. In other words, the core atoms 
begin to diffuse along axis and become 'wet', while the 
surface atoms remain 'solid-like' at this melting region. 
These results indicate that the melting mainly comes 
from the diffusion of core atoms and no surface melting 
happens at the beginning of the melting. Moreover, the 
dc fluctuates around 0.12 at low temperature, consistent 
with the Lindepiann criterion for equilibrium melting of 
simple crystalsEy, which shows that the interior melting 
takes place. For all the three cases, there is a rapid raise 
of rms bond length fluctuation in a narrow temperature 
region (1000-1150K), indicating that the surface atoms 
are involved into the melting process. Afterwards, the 
surface atoms play an important role on the melting of 
the wire at the high temperature region. Above 1150K, 
all of three quantities have a large, constant and smooth 
variation, corresponding to the completely melting sta- 
tus. Together with Fig. 2d, we hold the surface melting 
means the overall melting in the multi-walled ultrathin 
helical nanowires. The interior melting is in preference to 
the surface melting and the diffusion of core atoms has a 
dominant effect on the melting at low temperature. The 
surface atoms affect the melting at high temperature and 
surface melting represents the overall melting. 

The above interior melting behavior is obviously dif- 
ferent from that of the bulk, clusters and crystalline 
nanowires, where the surface melting usually occurs be- 
fore the total melting. In the bulk, after the melting 
of surface layers, the rest core atoms can still be seen 
as an analogous bulk. Much more energy is needed to 
make the 'rest bulk' molten. Therefore, the surface melt- 
ing temperature can be lower than the overall melting 



point. In the cases of clusters and crystalline nanowires, 
surface atoms have fewer nearest neighbors and are thus 
weaker bound and less constrained compared to the core 
atoms. Thus, surface atoms are easy to diffuse and be- 
come liquid-like. However, the present structures of gold 
nanowires are helical multi-walled cylindrical. The in- 
teraction among the same shell atoms is stronger than 
that of neighboring shell atoms. Furthermore, to reach a 
well-tightened helical multi-walled structure, the helical 
match may cause the fewer number of atoms in the inte- 
rior shells, especially for the center chain. Thus, interior 
atoms may have less coordinate numbers and larger inter- 
atomic distances than those of surface atoms. Therefore, 
core atoms in nanowires can break away from the binding 
sites in the wire prior to surface atoms and diffuse at a 
rather low temperature. In addition, in th optimization 
process of ground state structure, we have found that the 
formation of helical structure in the outer shell is earlier 
than the interior shell. This implies that the surface is 
dynamically more stable than the interior part. Exper- 
imentally, Wu et al. have found that the melting of Ge 
nanowire starts from, the two ends of the wire and move 
towards the middleE^. These imply that the melting be- 
havior for nanowiMjs dramatically differs from those for 
nanorods and bulkl3. 
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FIG. 3. The rms bond length fluctuation in the 
81(18-12-6-1) structure is ploted as functions of temperature 
for the entire system((5), the core atoms(5c) and the surface 
atoms((5s). 

To further clarify our ideas, we separate the function 
of surface atoms and core atoms by fixing the surface 
atoms artificially and allowing the interior atoms to move 
or vice verse. All the rms bond length fluctuation 5, Sfc, 
Sfs are calculated for the entire system, representing no 
constrained, fixing core atoms and fixing surface, respec- 
tively. We still take the wire SI as an example. As shown 
in Fig. 4, regardless the large difference among 6, Sfc, Sfs 
curves, the full melting temperatures (which corresponds 
to a large, smooth and constant rms bond length fluctu- 
ation) in the three cases are about 1150K. The Sfs has 
considerable jump at low temperature (400K), while the 
Sfc is almost invariable up to 800K. This indicates that 
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surface atoms with the hehcal structure are thermody- 
namically more stable. For the case of fixing the surface 
and no constrained, although the Sfs is just nearly half of 
the S, their general trends are similar. The small absolute 
value of Sfs comes from the fixing surface atoms, which 
contributes nearly zero to the rms bond length fluctua- 
tion. Therefore, we propose there should be no essential 
difference between the first two cases. These results also 
support the melting comes from the interior atoms at 
low temperature and the surface melting represents the 
overall melting. 
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FIG. 4. Melting behavior in the 18-12-6-1 structure with 
fixing interior atoms and no constraint. 

To check the validity of the current results, we have 
used the same MD code to study the melting of clus- 
ters and obtained xeaults similar to previous Monte Carlo 
(MC) simulationsLj'LJ. Wc have also exploited the melt- 
ing process of gold nanowires by MC and observed sim- 
ilar interior diffuse behavior. Moreover, to examine the 
effect of the periodic bound condition on the melting be- 
havior, we rescale the supercell length to 2-times and 
3-times of original one. The melting temperatures are 
also about HOOK and the observed melting process is 
similar to the above results. It proves that the periodic 
bound condition has little effect on the simulation re- 
sults. It should be further pointed out that we just limit 
our discussion to the helical structure in this paper. As 
mentioned above, this kind of helical structure is jpceya- 
lent in metal nanowires in the small diameter rangeli3~E3. 
Therefore, the current results on interior melting behav- 
ior in ultrathin gold nanowires is significant and might 
be a common feature in this kind of ultrathin metallic 
nanowires. 

In summary, the thermal behavior of helical multi- 
walled gold nanowires has been studied and the main 
points are made as following. (1) The melting process 
starts from the interior region and no surface melting 
happens at lower temperature. We further argued that 
interior melting behavior happens prior to surface melt- 
ing and surface melting represents the overall melting 
in ultrathin metallic nanowires. (2) The overall melting 



temperature of gold nanowires is lower than the one for 
the bulk, but higher than gold nanoclusters. (3) The 
surface and core atoms play different role in the melting 
behavior of these nanowires. The core atoms have a dom- 
inating effect on the melting at the beginning stage and 
the surface atoms are involved in the melting at higher 
temperature region. The core melting is closely related to 
the interior atomic structural characters. (4) The novel 
interior melting behavior is ultimately attributed to the 
helical multi- walled structure. 
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